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Recent, current, and planned NASA missions that employ blunt-body entry vehicles pose 
aerothermodyamic problems that challenge the state-of-the art of experimental and 
computational methods. The issues of boundary-layer transition and turbulent heating on 
the heat shield have become important in the designs of both the Mars Science Laboratory 
and Crew Exploration Vehicle. While considerable experience in these general areas exists, 
that experience is mainly derived from simple geometries; e.g. sharp-cones and flat-plates, or 
from lifting bodies such as the Space Shuttle Orbiter. For blunt-body vehicles, application 
of existing data, correlations, and comparisons is questionable because an all, or mostly, 
subsonic flow field is produced behind the bow shock, as compared to the supersonic (or 
even hypersonic) flow of other configurations. Because of the need for design and validation 
data for projects such as MSL and CEV, many new experimental studies have been 
conducted in the last decade to obtain detailed boundary-layer transition and turbulent 
heating data on this class of vehicle. In this paper, details of several of the test programs are 
reviewed. The laminar and turbulent data from these various test are shown to correlate in 
terms of edge-based Stanton and Reynolds number functions. Correlations are developed 
from the data for transition onset and turbulent heating augmentation as functions of 
momentum thickness Reynolds number. These correlation can be employed as engineering- 
level design and analysis tools. 
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Nomenclature 

skin friction coefficient 

vehicle diameter 

adiabatic wall enthalpy 

total free stream enthalpy 

wall static enthalpy 

wall static enthalpy at 300 K 

free stream Mach number 

boundary-layer edge Mach number 

free stream pressure 

Prandtl number 

heat transfer rate 

hemispherical nose radius 

free stream unit Reynolds Number PaJJ^/ 

free stream unit Reynolds Number p e U e / p e 

boundary-layer momentum thickness Reynolds Number p e U e G/p e 

Stanton number based on free stream conditions ^ w /[p oo f/ oo (^ 0 - H w )] 

Stanton number based on boundary layer edge conditions q w /[p e U e [H aw - H w )] 

free stream temperature 
free stream velocity 
boundary layer edge velocity 

normalized distance along model centerline for MSL local coordinate system 
normalized distance along model centerline for CEV local coordinate system 
angle of attack 

factor in laminar heating correlation [(Poo Poo )/( P e Pe {Pe/P™ )^ 4 
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turbulent heating augmentation factor SV measured /^predicted, laminar 
boundary layer momentum thickness 
free stream viscosity 
boundary layer edge viscosity 


I. Background 

Blunt body configurations are the most common geometries employed for entry into planetary atmospheres. 
Examples of manned blunt-body entry vehicles include the Mercury, Gemini, and Apollo capsules. Examples of 
unmanned flight test or interplanetary probe blunt-body entry vehicles are more numerous and include the Viking, 
Pioneer, FIRE II, ARD, OREX, Stardust, etc. An overview of major programs that includes vehicle and mission 
descriptions in which blunt body entry vehicles have been employed is given in Ref. 1 . 

Historically, turbulent flow at hypersonic speeds has generally been considered to be a design issue for moderate 
to high lift-to-drag ratio lifting body and winged vehicles, such as the Shuttle Orbiter, rather than for blunt body 
entry vehicles. These vehicles travel at high Reynolds numbers and their large size provides sufficient length over 
which turbulent flow can develop. In contrast, blunt-body re-entry vehicles generally are smaller and have high drag 
coefficients that cause rapid deceleration, which thus decreases the likelihood of producing turbulent flow. 

However, the problem of turbulent heating has become important in recent blunt-body vehicle designs, notably 
the Mars Science Laboratory (MSL) and the Orion Crew Exploration Vehicle (CEV). Both these vehicles are, in 
comparison to past blunt-body vehicles, quite large: the MSL is 4.5 m in diameter and the CEV will be ~5 m in 
diameter. Furthermore, both vehicles will experience atmospheric entry (at Mars and Earth, respectively) at high 
speeds (~1 1 km/s for CEV and ~5 km/s for MSL) and fly high angle-of-attack (for blunt bodies) lifting trajectories 
(1 1-deg to 16-deg range for MSL and 16-deg to 24-deg range for CEV). These factors all tend to promote transition 
to turbulent flow, and thus the conservative design philosophy applied in both projects is to assume turbulent flow 
throughout their trajectories. 

This assumption of turbulent flow led to the requirement to conduct high-Reynolds number hypersonic 
aeroheating testing on each vehicle in order to obtain turbulent heat transfer data. These data are used in the 
evaluation of the CFD models employed in the design of each vehicle. Several examples of these experimental 
studies will be presented in which comparisons with CFD results have been performed. 

II. MSL and CEV Project Descriptions 

In this discussion of blunt-body transition and turbulence, ground test programs and supporting computational 
analyses will be presented for the Mars Science Laboratory and the Project Orion Crew Exploration Vehicle. Brief 
overviews of each vehicle and mission are presented below. 

The Mars Science Laboratory (MSL) mission (Ref. 2), to be launched in 201 1, will deliver the largest (> 900 kg) 
rover ever to Mars. The MSL vehicle (Figure 1) will fly a controlled, lifting trajectory (a = 16 deg) to deliver the 
payload to within 10 km of the target location. The entry vehicle is comprised of a 4.5 m diameter spherically- 
blunted, 70-deg half-angle cone forebody heat shield that protects the aftbody payload from entry heating and 
provides a hypersonic lift-to-drag (L/D) ratio of 0.24 for aerodynamic control and maneuvering. The MSL vehicle 
will enter the Martian atmosphere at a velocity of 5.6 km/sec, which is greater than that of any other Mars probe 
except Mars Pathfinder. Because of the vehicle’s large heat shield diameter and the high angle of attack and 
velocity of its entry trajectory, it is expected that the vehicle will experience boundary-layer transition to turbulent 
flow well before the peak heating point on the trajectory. Therefore the vehicle’s Thermal Protection System (TPS) 
was designed (Refs. 3-4) with a conservative assumption of turbulent flow throughout the entire trajectory. 

The Project Orion Crew Exploration Vehicle (CEV) was defined by NASA’s Exploration Systems Architecture 
Study (Ref. 5) as NASA’s next manned space vehicle. The CEV will support NASA’s exploration missions by 
providing crew access to the International Space Station, the Moon, and Mars. The geometry of the CEV (Figure 2) 
is similar to that of Apollo - a spherical segment heat shield that protects a truncated-cone shaped crew 
compartment - but is considerably larger. The maximum diameter of the CEV is 5.0 m (current design), as 
compared to that of 3.912 m for Apollo. The design of the CEV’s TPS must account for the high heating rates 
generated at lunar return velocities and the aero thermo dynamic challenges of non-equilibrium thermo-chemistry, 
turbulent flow, and radiation transport. As with MSL, the CEV is being designed with the conservative assumption 
of fully-turbulent flow throughout its trajectory. Although not discussed herein, the CEV TPS design must also 
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account for ablation effects, including shape change, flow field radiation-ablation coupling, and ablated surface 
roughness heating augmentation. An overview of the CEV design and development program is given in Ref 6 and a 
summary of all CEV aerothermodynamic testing to date is presented in Ref. 7. 
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Figure 1. Mars Science Laboratory Entry Vehicle 



III. Survey of Blunt Body Turbulent Experimental Data Sources 

Extensive aerothermodynamic ground-testing has been conducted in support of both the MSL and CEV 
development program developments, including testing in the AEDC Hypervelocity Tunnel 9, CUBRC LENS I and 
48-Inch Reflected Shock Tunnels, CalTech T5 Reflected Shock Tunnel, and the NASA Langley 20-Inch Mach 6 Air 
and 31 -Inch Mach 10 Air Tunnels. In the analyses presented herein, only the test data obtained on smooth outer- 
model line (no penetrations, protrusions, roughness, or other effects, except for boundary-layer trips) models at 
perfect-gas conditions where transitional and/or turbulent data flow was produced were considered. These tests are 
summarized in Table 1. Nominal flow conditions for the data sets under consideration are listed in Table 2 for the 
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AEDC Hypervelocity Tunnel 9, LENS I, and the NASA Langley 20-Inch Mach 6 Air Tunnel. Information on these 
facilities can be found in Refs. 8-10. 

A brief discussion of each test entry along with sample data and comparisons will be presented below. More 
detailed results on model design, run schedules and exact test conditions, along with the full data sets, can be found 
in the references given. In subsequent sections, analyses and correlations for transition onset and turbulent heating 
augmentation that are based on the entire data set from all of the tests will be presented. 


Table 1. Summary of CEV and MSL Tests with Transitional / Turbulent Data 


Test 

Facility 

Test 

Number 

Config. 

Model 

diam. 

a (deg) 

Model 

fabrication 

Instrumentation 

Refs. 

Notes 

AEDC Tunnel 
9, Mach 8 and 
Mach 10 

N/A 

MSL 

6-in. 

0-24 

Stainless 

steel 

Coaxial surface 
thermocouples 

11 

High Reynolds 
number turbulent 
aeroheating 

LaRC 20-Inch 
Mach 6 Air 

Test 6827 

MSL 

5-in. 

11, 16, 
20 

Cast 

ceramic 

Global phosphor 
thermography 

12, 13, 14 

Boundary-layer 

tripping 

LaRC 20-Inch 
Mach 6 Air 

Test 6884 

MSL 

6-in. and 
7-in. 

16, 20 

Cast 

ceramic 

Global phosphor 
thermography 

15 

Transition and 
turbulence 

LaRC 20-Inch 
Mach 6 Air 

Test 6945 

MSL 

7-in., 8- 
in. and 9- 
in. 

16 

Cast 

ceramic 

Global phosphor 
thermography 

Herein 

Transition and 
turbulence 

AEDC Tunnel 
9, Mach 8 and 
Mach 10 

39-CH 

CEV 

7-in. 

24-32 

Stainless 

steel 

Coaxial surface 
thermocouples 

16, 17 

High Reynolds 
number turbulent 
aeroheating 

CUBRC 
LENS I 

67-CH 

CEV 

14-in. 

20, 28 

Stainless 

steel 

Coaxial surface 
thermocouples and thin- 
film heat transfer gages 

18 

High Reynolds 
number turbulent 
aeroheating 

LaRC 20-Inch 
Mach 6 Air 

Test 6917 
(31-CH) 

CEV 

5 -in., and 
7-in. 

28 

Cast 

ceramic 

Global phosphor 
thermography 

19 

Boundary-layer 

tripping 

LaRC 20-Inch 
Mach 6 Air 

Test 6931 
(56-CH) 

CEV 

7-in. 

16-32 

Stainless 

steel 

Coaxial surface 
thermocouples 

20,21 

Boundary-layer 

tripping 

LaRC 20-Inch 
Mach 6 Air 

Test 6944 

CEV 

7-in., 8- 
in. and 9- 
in. 

28 

Cast 

ceramic 

Global phosphor 
thermography 

herein 

Transition and 
turbulence 


Table 2. Facility Test Conditions 


Test 

Facility 

Test 

Gas 

Re. 

(1/ft) 

Re. 

(1/m) 

M. 

p. 

(Pa) 

T. 

(K) 

p, 

(kg/m A 3) 

U. 

(m/s) 

H 0 -H 300 K 

(MJ/kg) 

AEDC Tunnel 9, Mach 8 

N2 

4.09E+06 

1.34E+07 

7.4 

1171.1 

74.4 

5.305E-02 

1302.2 

0.536E+06 

AEDC Tunnel 9, Mach 8 

N2 

8.35E+06 

2.74E+07 

7.4 

2341.2 

73.5 

1.073E-01 

1301.1 

0.591E+06 

AEDC Tunnel 9, Mach 8 

N2 

16.4E+06 

5.37E+07 

7.6 

5177.9 

81.1 

2.152E-01 

1399.1 

0.692E+06 

AEDC Tunnel 9, Mach 8 

N2 

21.6E+06 

7.08E+07 

7.8 

6207.3 

77.1 

2.716E-01 

1390.0 

0.674E+06 

AEDC Tunnel 9, Mach 8 

N2 

30.5E+06 

10.0E+07 

7.8 

8236.8 

73.8 

3.766E-01 

1360.7 

0.671E+06 

AEDC Tunnel 9, Mach 8 

N2 

48.6E+06 

16.0E+07 

8.0 

11855.3 

70.1 

5.718E-01 

1356.8 

0.681E+06 

AEDC Tunnel 9, Mach 10 

N2 

1.22E+06 

0.40E+07 

9.6 

158.2 

51.8 

1.029E-02 

1401.7 

0.671E+06 

AEDC Tunnel 9, Mach 10 

N2 

1.90E+06 

0.62E+07 

9.6 

272.9 

55.6 

1.653E-02 

1459.8 

0.781E+06 

AEDC Tunnel 9, Mach 10 

N2 

4.67E+06 

1.53E+07 

9.9 

606.7 

52.9 

3.863E-02 

1463.0 

0.803E+06 

AEDC Tunnel 9, Mach 10 

N2 

8.79E+06 

2.88E+07 

10.1 

1092.8 

52.1 

7.072E-02 

1482.1 

0.822E+06 

AEDC Tunnel 9, Mach 10 

N2 

14.2E+06 

4.66E+07 

10.3 

1588.1 

49.3 

1.086E-01 

1475.0 

0.818E+06 

AEDC Tunnel 9, Mach 10 

N2 

20.0E+06 

6.56E+07 

10.4 

2175.2 

48.5 

1.51 IE-01 

1472.2 

0.805E+06 

LaRC 20-Inch Mach 6 Air 

Air 

2.13E+06 

0.70E+07 

6.0 

586.7 

61.9 

3.301E-02 

939.1 

0.202E+06 

LaRC 20-Inch Mach 6 Air 

Air 

2.61E+06 

0.86E+07 

6.0 

723.2 

62.2 

4.050E-02 

943.0 

0.206E+06 

LaRC 20-Inch Mach 6 Air 

Air 

2.99E+06 

0.98E+07 

6.0 

825.0 

62.2 

4.624E-02 

945.5 

0.209E+06 

LaRC 20-Inch Mach 6 Air 

Air 

3.44E+06 

1.13E+07 

6.0 

935.5 

61.6 

5.290E-02 

940.1 

0.202E+06 

LaRC 20-Inch Mach 6 Air 

Air 

3.88E+06 

1.27E+07 

6.0 

1098.0 

63.3 

6.052E-02 

957.2 

0.221E+06 

LaRC 20-Inch Mach 6 Air 

Air 

4.13E+06 

1.36E+07 

6.0 

1118.5 

61.5 

6.345E-02 

942.1 

0.205E+06 

LaRC 20-Inch Mach 6 Air 

Air 

4.63E+06 

1.52E+07 

6.0 

1321.7 

63.5 

7.250E-02 

985.5 

0.248E+06 

LaRC 20-Inch Mach 6 Air 

Air 

5.09E+06 

1.67E+07 

6.0 

1441.5 

63.4 

7.920E-02 

985.7 

0.248E+06 

LaRC 20-Inch Mach 6 Air 

Air 

5.76E+06 

1.89E+07 

6.0 

6474.7 

63.0 

3.590E-01 

957.1 

0.220E+06 
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Test 

Facility 

Test 

Gas 

Re. 

(1/ft) 

Re. 

(1/m) 

M„ 

p. 

(Pa) 

T. 

(K) 

p. 

(kg/m A 3) 

Uoo 

(m/s) 

HcrPfiociK 

(MJ/kg) 

LaRC 20-Inch Mach 6 Air 

Air 

6.63E+06 

2.18E+07 

6.0 

1833.0 

62.6 

1.023E-01 

954.6 

0.217E+06 

LaRC 20-Inch Mach 6 Air 

Air 

7.32E+06 

2.40E+07 

6.1 

2032.0 

62.7 

1.130E-01 

956.8 

0.220E+06 

LaRC 20-Inch Mach 6 Air 

Air 

7.46E+06 

2.45E+07 

6.0 

2058.0 

62.5 

1.154E-01 

953.6 

0.217E+06 

CUBRC LENS I 

Air 

8.65E+05 

0.28E+07 

7.6 

260.8 

79.8 

1.134E-02 

1371.4 

0.719E+06 

CUBRC LENS I 

Air 

9.40E+05 

0.31E+07 

7.7 

231.0 

70.6 

1.134E-02 

1304.5 

0.620E+06 

CUBRC LENS I 

Air 

9.51E+05 

0.31E+07 

7.7 

238.3 

71.5 

1.154E-02 

1312.7 

0.632E+06 

CUBRC LENS I 

Air 

7.92E+06 

2.60E+07 

7.7 

2178.7 

75.8 

9.947E-02 

1353.3 

0.691E+06 

CUBRC LENS I 

Air 

23.4E+06 

7.68E+07 

8.0 

6095.0 

78.9 

3.036E-01 

1420.4 

0.787E+06 

CUBRC LENS I 

Air 

25.8E+06 

8.46E+07 

8.1 

5364.1 

65.6 

2.835E-01 

1316.7 

0.631E+06 

CUBRC LENS I 

Air 

26.3E+06 

8.63E+07 

8.1 

5968.6 

69.1 

2.991E-01 

1346.2 

0.674E+06 

CUBRC LENS I 

Air 

30.4E+06 

9.96E+07 

8.2 

7411.9 

73.1 

3.525E-01 

1400.6 

0.753E+06 

CUBRC LENS I 

Air 

32.5E+06 

10.7E+07 

8.0 

8756.3 

76.1 

3.979E-01 

1396.0 

0.749E+06 


IV. MSL Turbulent Aeroheating Data 

A. MSL Turbulent Aeroheating in AEDC Hypervelocity Tunnel 9 

An investigation (Ref. 11) of turbulent aeroheating on the MSL vehicle was conducted in the Arnold 
Engineering Development Center (AEDC) Hypervelocity Wind Tunnel No. 9. In this study, aeroheating data were 
collected on a 6-in. (0.1524 m) diameter coaxial thermocouple instrumented MSL model in perfect-gas N 2 flow at 
the tunnels Mach 8 and Mach 10 test conditions. Data were obtained at free stream Reynolds numbers of 4xl0 5 6 /ft to 
49xl0 6 /ft at Mach 8 and at lxl0 6 /ft to 19xl0 6 /ft at Mach 10 with angles-of-attack between 0-deg and 24-deg. 
Turbulent flow was produced over the leeside of the heat shield at the highest Mach 10 Reynolds number and 
turbulent flow was produced over the entire heat shield (both leeside and wind-side) at the highest Mach 8 Reynolds 
number. The experimental uncertainty of the data was estimated to be ±12%. 

Laminar and turbulent perfect-gas comparisons to these data were performed using the LAURA code (Ref. 22) 
with the algebraic Cebeci-Smith turbulent model for the turbulent cases. Comparisons between centerline data and 
predictions are shown for selected a = 16-deg cases in Figure 3 and Figure 4. For the Mach 10 cases, the laminar 
predictions and data were in close agreement at Re.= 4xl0 6 /ft case, while at Re = 19xl0 6 /ft, the laminar predictions 
matched the data on the wind-side of the forebody and the turbulent predictions matched the data on the leeside. For 
both Mach 8 cases, the turbulent predictions matched the data. Although the agreement between predictions and 
data was generally good for these cases, there were discrepancies around the stagnation point ( x/R ~ -0.4) for all 
cases where measured heating rates were considerably higher than either laminar or turbulent predictions. 



Figure 3. MSL Data and Comparisons from AEDC Tunnel 9 at Mach 10, a=16-deg 
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Figure 4. MSL Data and Comparisons from AEDC Tunnel 9 at Mach 8, a=16-deg 


B. MSL Transitional/Turbulent Aeroheating in the NASA Langley 20-Inch Mach 6 Air Tunnel 

Several test programs were performed in the Langley 20-Inch Mach 6 Air Tunnel. In all cases, testing was 
performed with global thermography on phosphor-coated ceramic MSL models. The experimental uncertainty of 
the resulting data was estimated to be ±10-15%. 

7. Test 6827 - Cavity effects on heating and transition 

An extensive test program (Refs. 12, 13, 14) was conducted during the configuration development phase of the 
MSL program to determine the heating augmentation and boundary- layer transition influence of circular cavities on 
the vehicle’s heat shield. These cavities were attachment points to the spacecraft cruise module that were located on 
the heat shield in early version of the MSL design, but were subsequently moved to the aftbody of the entry vehicle 
where they will have minimal effects on heating. Data were obtained at free stream Reynolds numbers of 2.0x1 0 6 /ft 
to 5.8xl0 6 /ft at angles of attack of 11, 16, and 20-deg on 5-in. diameter models. Although the majority of these tests 
were conducted on models with the attachment point cavities, a limited number of runs were made using smooth 
surface models that had an array of boundary-layer trips located to the leeside of the nose. These previously un- 
reported runs were part of Test 6827 and are shown here for the first time for use in the current analysis. The 
transitional / turbulent data from these runs are plotted in Figure 5 along with comparison to laminar LAURA 
predictions (no turbulent predictions were performed). The trips appeared to produce fully -turbulent flow when 
located at the nose/cone tangency point but not when located midway down the conical section. 

2. Test 6884 - Transition and turbulence with cone-angle parametrics 

In Test 6884, a parametric investigation of the effects of cone angle on transition and turbulent heating was 
conducted. MSL 70-deg sphere-cone models of 6-in. and 7-in. (along with 6-in. diam. 50-deg and 60-deg. sphere- 
cone models) were tested at free stream Reynolds numbers of 2.0xl0 6 /ft to 7.3xl0 6 /ft at angles of attack of 11, 16 
and 20-deg. Sample results are shown in Figure 6 and addition data are presented in Ref. 15. Transitional flow 
occurred on both size models, but fully-developed turbulent flow was not produced. 

3. Test 6945 - Model size effects on blockage and transition 

The inability to generate turbulent flow on the MSL configuration without using boundary-layer trips in previous 
tests led to an investigation (Test 6945) into the possibility of using larger diameter models to produce turbulent 
flow without adversely affecting the test section flow field. MSL models of 7-in., 8-in. and 9-in. model diameters 
were tested at a =16-deg at free stream Reynolds numbers of 2.0x1 0 6 /ft to 7.5xl0 6 /ft. Data from this test are 
presented for the first time here in Figure 7 and comparisons with LAURA predictions are shown in Figure 8. It was 
found the 9-inch model caused definite flow blockage and thus good quality data could not be obtained on this 
model, but the 8-inch model appeared to have no effect on flow quality and turbulent flow was observed on the outer 
leeside portion of the cone. However, the data were approximately 15% to 20% higher than turbulent predictions. 
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Stx(Re soD ) (1/2) Stx(Re, D ) 


Possibly these difference were due to the fact that the solutions were generated with turbulent flow over the entire 
body, whereas the data indicate that transition occurred just downstream of the nose/cone tangency point. The 7-in. 
diameter model data were similar to those from previous tests. 
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LAURA: Laminar 

V Test 6827, Run 76 a=16, Re,=2.6x10 6 /ft 

► Test 6827, Run 75 a=16, Re,=3.1x10 6 /ft 

< Test 6827, Run 74 a=16, Re,=3.4x10 8 /ft 

♦ Test 6827, Run 73 «=16, Re, =4.2x1 0 6 /ft 

• Test 6827, Run 72 a=16, Re,=5.1x10 6 /ft 

■ Test 6827, Run 70 u=16, Re,=5.8x10 6 /ft 
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LAURA: Laminar 

► Test 6827, Run65a=16, Re,=3.1 x10 6 /ft 
< Test 6827, Run 66 a=1 6, Rel=3.4x1 0 6 /ft 
♦ Test 6827, Run 67 a=16, Re x =4.2x10 6 /ft 
O Test 6827, Run68a=16, Re^=5.1 x10 6 /ft 
■ Test 6827, Run 69 a=1 6, Re"=5.8x1 0 6 /ft 



5-inch diameter model, a =16-deg, trip array at 
nose/cone tangency point 


5 -inch diameter model, a =16-deg, trips midway along 
conical section 


Figure 5. MSL cx=16-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6827 




6-inch diameter model, a =16-deg 


x/R 

7-inch diameter model, a =16-deg 


Figure 6. MSL a=16-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6884 
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Figure 7. MSL Global Heating Data from Test 6945: 8-inch diameter model 
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Figure 8. MSL a=16-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6945 
V. CEV Turbulent Aeroheating Data 

A. CEV Turbulent Aeroheating at AEDC Hypervelocity Tunnel 9 

Turbulent aeroheating testing (Refs. 16, 17) of the Orion CEV vehicle was conducted in the AEDC 
Hypervelocity Wind Tunnel No. 9 on a 7-inch diameter, thermocouple-instrumented model in perfect-gas N 2 flow at 
the tunnels Mach 8 and Mach 10 test conditions. Angles-of-attack were varied from 16-deg to 32-deg. Free stream 
Reynolds numbers ranged from 8xl0 6 /ft to 48xl0 6 /ft at Mach 8 and from 2xl0 6 /ft to 20xl0 6 /ft at Mach 10. 
Transitional or turbulent flow was produced on the heat shield leeside for Re* > 10xl0 6 /ft for both Mach 8 and Mach 
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10, and for the higher Mach 8 Reynolds numbers, turbulent flow was also produced on the wind side of the heat 
shield. The experimental uncertainty of the data was estimated to be ±12%. 

Laminar and turbulent perfect-gas comparisons to these data were performed using the LAURA code with the 
algebraic Cebeci- Smith turbulent model being used for the turbulent cases. Comparisons between centerline data 
and predictions are shown for selected a = 28-deg cases in Figure 9 and Figure 10. For the Mach 10 cases, the 
laminar predictions and data were in close agreement at the R e. = 4xl0 6 /ft case, while at Re*= 19xl0 6 /ft, the laminar 
predictions matched the data on the wind-side of the forebody and the turbulent predictions matched the data on the 
leeside. For both Mach 8 cases, the turbulent predictions matched the data. 



Figure 9. CEV Mach 10, a=28-deg Data and Comparisons from AEDC Tunnel 9 



Figure 10. CEV Mach 10, a=28-deg Data and Comparisons from AEDC Tunnel 9 
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B. CEV Turbulent Aeroheating at CUBRC LENS 

Additional high Reynolds number turbulent aeroheating testing of the CEV was performed in the CUBRC LENS 
I reflected shock tunnel (Ref. 18) to supplement the AEDC test with data at lower (a =20-deg) angle of attack. Data 
were obtained with a mix of coaxial surface thermocouples and thin- film heat- transfer gages on a 14-in. diam. 
stainless steel model. Runs were performed at free stream Reynolds numbers of 0.9x1 0 6 /ft to 32xl0 6 /ft. Fully- 
turbulent flow was produced over the entire forebody at the highest Reynolds numbers and close agreement between 
predictions and data were obtained, as shown in Figure 11. 



Mach 8, Re* = 0.96x1 0 6 /ft, laminar data Mach 8, Re* = 32xl0 6 /ft, turbulent data 


Figure 11. CEV Mach 8, a=20-deg Data and Comparisons from CUBRC LENS Test 67-CH 


C. CEV Transitional/Turbulent Aeroheating in the NASA Langley 20-Inch Mach 6 Air Tunnel 

Several CEV tests were performed in the Langley 20-Inch Mach 6 Air Tunnel using both stainless-steel, 
thermocouple-instrumented models and phosphor-coated, cast ceramic models with global thermography. The 
experimental uncertainty of the resulting data was estimated to be ± 10% to 15%. 

7. Test 691 7 - Boundary -layer tripping 

In Test 6917, turbulent flow was produced on the leeside of the 7-in. diam. CEV models using arrays of 
boundary layer trips. Data were obtained at a =28-deg at free stream Reynolds numbers of 2.1xl0 6 /ft to 7.3x1 0 6 /ft 
as reported in Ref. 19. Laminar and turbulent (with the Cebeci-Smith algebraic model) predictions were performed 
using LAURA and close agreement was obtained as shown in Figure 12. 

2. Tests 6931 and 6932 - Boundary -layer tripping 

In Test 6931 and 6932, the coaxial-instrumented, 7-in. diam. CEV model tested at AEDC Tunnel 9 was tested in 
the Langley 20-Inch Mach 6 Air Tunnel. Test 6931 was a checkout of the model (Ref. 20) and Test 6932 (Ref. 21) 
was a study of the effects of boundary-layer tripping on both the forebody and aftbody of the model. Runs were 
performed in both static and continuous pitch-mode over an angle-of-attack range of 16-deg to 32-deg free stream 
Reynolds numbers of 2.1xl0 6 /ft to 7.3x1 0 6 /ft. Sample data and comparisons are given in Figure 13. 

3. Test 6944 - Model size effects on blockage and transition 

In conjunction with the model size blockage and transition testing conducted for MSL (Test 6945), a similar 
investigation was performed for the CEV configuration. CEV models of 7-in., 8-in. and 9-in. diameters were tested 
at a =28-deg at free stream Reynolds numbers of 2.0xl0 6 /ft to 7.5xl0 6 /ft. These Test 6944 data are presented for 
the first time herein in Figure 14. As with the MSL blockage study, the 9-in. diam. CEV model caused tunnel 
blockage and unsteady flow conditions. However, the 7-in. and 8-in. diameter models were successfully run and 
fully turbulent flow was produced on the leeside of the forebody at the highest test Reynolds numbers as shown in 
Figure 15. 


10 

American Institute of Aeronautics and Astronautics 



z/R 



z/R 


7-inch diameter model, a = 28-deg, Re. = 5.6x 10 6 /ft 7-inch diameter model, a = 28-deg, Re. = 7.4x 10 6 /ft 


Figure 12. CEV a=28-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6917 




7-inch diameter model, a = 28-deg, trip array on 7-inch diameter model, a = 28-deg, trip array on negative 
positive side of symmetry axis side of symmetry axis 


Figure 13. CEV a=28-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6932 
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Re = 3.0E6/ft 


Re = 3.9E6/ft 


Re = 5.0E6/ft 



Figure 14. CEV Global Heating Distribution from Test 6944, a = 28-deg, 8-in. diam. model 
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Mach 6, 8-inch diameter model, a = 28-deg 


Figure 15. CEV a=28-deg Data and Comparisons from LaRC 20-Inch Mach 6 Air Test 6944 

VI. Comparison and Correlation of Heating Data Sets 

Having presented an overview of the various experimental turbulent aeroheating data sets available on blunt 
bodies in the previous section, facility-to-facility correlations of the data will now be demonstrated. Such 
correlations can provide a valuable check on experimental methods employed in each test and the resulting data. 


12 

American Institute of Aeronautics and Astronautics 


In the heating plots shown previously herein, the data have been presented in terms of a correlation parameter 
based on free stream conditions: 


( 1 ) Stx^Re 0 


,d 


pJJ„{H 0 -H w ) 


pooU^D 


Pa 


As has been shown in these plots, this parameter correlates the Reynolds number dependence of the laminar 
heating distributions to a separate, nearly-constant curve for each facility at a given angle-of-attack. However, close 
examination of these figure reveals that the correlation is facility-dependent, as shown in Figure 16. The reason that 
these correlations are not universal is that the laminar heating rates have Mach number compressibility and enthalpy 
ratio dependencies and these parameters vary from facility to facility. However, a universal form of the correlation 
can be determined by reference to the original boundary- layer based heat-transfer correlation form. 

The St x (Re. ?D )° 5 factor can be derived from boundary- layer correlations for stagnation point heating, such as 
those by Fay and Riddell or Van Driest (e. g. Ref. 23). Since the data in question are all from low-enthalpy, perfect- 
gas test conditions, the Van Driest formulation can be employed for simplicity: 
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With some algebraic re-arrangement and the use of the Euler / Newtonian approximations for the velocity 
gradient term, Eq. (1) becomes 
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where the Stanton and Reynolds numbers are defined in terms of boundary-layer edge conditions, not free 
stream conditions. With further algebra, this equation can be re-formulated in terms of free stream conditions as: 
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For perfect-gas, stagnation point conditions, the assumptions can be made that Pr = constant and H aw = H 0 . Thus: 
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Now define a laminar correlation factor as 
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For a hemisphere, R n = D/2, but for spherical arc-segment (e.g. CEV) or sphere-cone (e.g. MSL) an effective nose 
radius is employed in place of a hemispherical nose radius, where: 

(7) Ru,eff = 7) x C shape _ factor ; where C shape _ factor is a geometry - and angle of attack - dependent constant 
Then finally 

(8) ^St x ^Re^ D j = local constant = f(geometry,a) 

It should also be noted that the original 0.763 leading constant in Eq. (2) is actually dependent on the 
atmospheric compositions (e.q. Ref. 24). However, for the perfect-gas Air and N 2 of these data sets, that constant 
differs by less than 1%. 

While this analysis is derived from a stagnation-point equation, the flow over the entire surface of a blunt body is 
similar enough to stagnation point flow that Eq. (8) applies anywhere on the body, albeit with different values for 
the local constant. That local constant could, if explicitly desired, be determined from the wind tunnel data (or 
estimated by various theoretical methods), but the purpose of this analysis was to determine a formulation by which 
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the various data sets could be correlated. And, in contrast to the free stream parameter correlations in Figure 16, the 
laminar data shown in Figure 17 are well correlated in terms of the edge-condition based definition of Eq. (8). The 
heating distributions from a wide range of free stream Reynolds numbers in three different facilities (or four 
facilities if the AEDC Tunnel 9 Mach 8 and Mach 10 nozzles are considered separately) have been collapsed to a 
single band of data with a spread of about ±15% (except around the stagnation point, where difference from laminar 
predictions have been noted in several studies - e.g. Refs. 10, 17). 

It should be noted that although these edge-based parameters provide a much better correlation for the MSL and 
CEV configurations, they are typically not used in standard data-reduction schemes. The free-stream based 
parameters are simpler to determine and do tend to correlate the data from a single facility. Additionally, the edge- 
based definitions contain assumptions that apply only to a blunt-body stagnation region and would not be applicable 
to a flat-plate, small-angle sphere-cone, or complex geometry - see for example the correlation for small-angle 
biconics in Ref. 25 which includes free-stream Mach number and cone-angle corrections. Thus, to avoid ambiguity, 
data from wind tunnel are typically defined in terms of free stream conditions such as shown in Figure 3 - Figure 15. 


-■ AEDC Tunnel 9, Mach 10 

-▼ AEDC Tunnel 9, Mach 8 

LaRC 20-Inch Mach 6, Test 6945 


CEV a=28-deg data 


MSL a=16-deg data 




Figure 16. Failure of Free Stream Parameters to Correlate Laminar Data from Different Facilities 
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Figure 17. Correlation of Laminar Data Using Corrected Parameters 
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While the edge-based definitions can be used to correlate laminar data, they fail when applied to the turbulent 
data as is shown in Figure 18. 
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Figure 18. Failure of Corrected Laminar Parameters to Correlate Turbulent Data 

This result is to be expected as the original stagnation-point correlation does not include the effects of 
turbulence. Employing an analogy to incompressible boundary-layer flow over a flat-plate, the laminar and 
turbulent skin friction coefficients (and by Reynolds analogy, the heat-transfer coefficients) have the forms: 
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Based on this relations, it could be hypothesized that a turbulent correlation function analogous to the laminar 
correlation function would be given by: 
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Unfortunately, evaluation of the data sets in terms of these parameters did not result in a valid correlation of the 
turbulent data. Instead, through trial-and-error empiricism, it was found that the turbulent data could be correlated 
with slightly different exponents: 
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This formulation was applied to all the turbulent CEV and MSL from all tests. As shown in Figure 19 and 
Figure 20, the turbulent data for each angle-of-attack correlated to within ~ ± 15%, including both the natural and 
tripped transition cases. Correlation formulations have thus been validated for data from multiple facilities and test 
conditions for laminar heating - as given by Eq. (8) - and for turbulent heating - as given by Eq. (12). 
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CEV a=20-deg data 
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Figure 19. Correlation of CEV Turbulent Heating Data Using Corrected Turbulent Parameters 
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Figure 20. Correlation of MSL Turbulent Heating Data Using Corrected Turbulent Parameters 


VII. Turbulent Heating Augmentation and Transition Onset Correlations 

For configuration analysis and design trade studies, it is useful to have engineering-level approximation for the 
onset of boundary-layer transition and the for resulting turbulent augmentation above laminar levels. While 
transition is vastly-complex process, engineering literature provides numerous examples of configuration and 
condition dependent transition correlations of reasonable accuracy. In that spirit, correlations were sought for the 
MSL and CEV configurations with the limiting assumptions that bound the data sets considered: low-enthalpy, 
perfect-gas air or N 2 flows; cold-wall conditions; subsonic boundary- layer edge conditions; minimal surface 
roughness; no blowing or ablation; and minimal cross-flow effects. 

The response function considered in the generation of these correlations was the ratio of the measured Stanton 
number to the predicted, cold- wall laminar value. 

C f 

(14) ^ measured 

^predicted , laminar 
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A predicted value for the dominator was chosen to provide a consistent definition and avoid any statistical 
fluctuations in the measurements that could bias the results. Based on the comparisons shown herein and the listed 
references, it was found that the agreement between predicted and measured laminar heating levels for the CEV and 
MSL geometries was within less than ± 10% and thus the predictions were suitable for use in correlations 

With respect to the experimental data employed in the numerator, it is emphasized that the domain of the data set 
includes all angle-of-attack and all Reynolds numbers from all facilities, the only exception being the tripped 
boundary-layer data. Furthermore the heating data were considered for the entire forebody (not just along the 
centerline) up to the rapid flow-expansion point ahead of the vehicle shoulder. A very large data set was thus 
considered since the AEDC thermocouple models included 38 off-centerline gages for the CEV geometry and 12 
off-centerline gages for the MSL geometry and the global thermographic phosphor data sets from the LaRC 20-Inch 
Mach 6 Air tunnel tests contain 512 x 480 pixel arrays. 

Following commonly employed practice for transition and heating correlations, the boundary-layer parameters 
Re 0 and Re 0 /M e were considered as forcing functions. These values were extracted from the predicted flow fields 
using a boundary- layer edge location definition as the point where the local enthalpy is 0.995 of the free stream total 
enthalpy. 

The resulting correlations for MSL and CEV in terms of Re 0 /M e are presented in Figure 21. Given that these 
flow fields all have subsonic to barely supersonic (M e < 1.5) edge conditions, this method was not expected to 
produce good results. And, in fact, the CEV data were completely uncorrelated in terms of Re 0 /M e . The MSL 
results did though provide a reasonable correlation: transition onset occurs between 300 < Re 0 /M e < 600 and the 
augmentation factor asymptotes to a linear relationship with Re 0 /M e . 

As shown in Figure 22, the correlations in terms of Re e vs. § were considerably better than those in terms of 
Re 0 /M e vs ((). Transition onset for the CEV geometry is clearly defined at around Re 0 = 200, although a slight inverse 
relationship between free stream Reynolds number and transition Re 0 was observed. Whether this is a truly flight- 
relevant relationship or a function of increasing free stream noise levels with Reynolds number in the wind tunnel 
testing is unclear. For the MSL geometry, a slightly less defined Re 0 transition onset was observed (although 
considerably better defined that for Re 0 /M) with values between Re 0 = 200 to 300. A similar inverse relationship was 
observed between free stream Reynolds number and transition onset Re 0 . 

The complete failure of the Re 0 /M e formulation for the CEV dataset but not for the MSL dataset can be explained 
by the relative complexity of the surface geometries and resulting flow fields. The CEV is a simple spherical-cap 
segment, whereas the MSL is a sphere-cone. Samples of the resulting Re 0 /M e and Re 0 distributions produced over 
these different geometries are shown in Figure 23 and Figure 24. While Re 0 varies over the surface for both 
geometries, Re 0 /M e is very nearly constant over the CEV geometry; thus it cannot be used as a correlation parameter 
for the CEV. 
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Figure 21. Correlation of Turbulent Heating Augmentation Factor in Terms of Re 0 /M e 
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CEV Data MSL Data 

Figure 22. Correlation of Turbulent Heating Augmentation Factor in Terms of Re 0 



CEV, a=28-deg MSL, cx=16-deg 

Figure 23. Sample Re 0 distributions for CEV and MSL, Mach 6, Re.=7.5xl0 6 /ft 
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Figure 24. Sample Re 0 /M e distributions for CEV and MSL, Mach 6, Re„=7.5xl0 6 /ft 
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With respect to transition, these results tend to confirm commonly-employed “rule-of-thumb” assumptions that 
transition on smooth blunt bodies occurs for an Re e of approximately 200. With respect to heating, the correlations 
show that for both CEV and MSL, the turbulent heating augmentation factors asymptote to linear functions of Re 0 
and 0 given by: 

(15) 0 = 1 + — J (for CEV) 
and 

R P 

(16) 0 = 1 + — (for MSL) 

158 

These two functions are similar enough that to suggest that for subsonic flow, there may be a universal 
correlation function for blunt-body turbulent heating augmentation. However, more data would be required on other 
blunt-body geometries to validate such an assumption. 

VIII. Summary and Recommendations 

An overview of recent turbulent aeroheating studies on the CEV and MSL blunt body geometries has been 
presented. Sample data and comparisons have been presented for each test. From these data, edge-based parameters 
have been validated for the correlation of both laminar and turbulent heating data from multiple facilities and test 
conditions. These parameters are: 
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Transition onset and turbulent heating augmentation factors were also determined from these data and correlated 
with the boundary-layer momentum thickness Reynolds number, Re 0 . For the CEV and MSL geometries, it was 
found that transition onset was given by: 

(19) Re 0 « 200 to 300 

and the turbulent heating augmentation factors asymptotically approached limits given by: 

(20) 0-1 + ^ (for CEV) 

R P 

(21) 0 = 1 + (for MSL) 

1 JO 

The relationships for laminar and turbulent heating and boundary- layer transition onset given by Eqs. (17) - (21) 
provide engineering-level design methods that can be employed in configuration design, development and 
optimization studies of blunt-body entry vehicles. 

These functions have been shown to be provide good representations of the data from which they were derived; 
however, significant gaps in the extent of the data can be inferred from consideration of the test conditions. 
Specifically, the experimental data sets do not encompass the effects on transition and heating of: high-enthalpy 
chemistry (total enthalpy levels » 1 MJ/kg); “hot” walls (T w » 300 K); test-gas effects (e.g. C0 2 for Mars); or 
surface roughness and blowing (due to ablation). These effects have not been considered because, for blunt bodies 
such as CEV and MSL, there are very few experimental data sets available in which these phenomena are present. It 
is thus recommended that acquisition of such data should be a priority in future research. 
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